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Abstract. We consider the inverse problem of the broken ray transform (sometimes 
^ , also referred to as the V-line transform). Explicit image reconstruction formulas are 

derived and tested numerically. The obtained formulas are generalizations of the 
filtered backprojection formula of the conventional Radon transform. The advantages 
of the broken ray transform include the possibility to reconstruct the absorption and 
t"*"- ' the scattering coefficients of the medium simultaneously and the possibility to utilize 

^^ , scattered radiation which, in the case of the conventional X-ray tomography, is typically 

discarded. 
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1. Introduction 

Image reconstruction tecliniques based on inversion of the Radon transform are 
well-established. While the classical Radon transform utilizes straight rays, there 
exists a considerable interest in the physical situations wherein the straight-ray 
propagation is lost due to the effects of scattering or refraction, yet measurements 
can be mathematically related to integrals of the medium properties over well-defined 
trajectories. This leads to various generalizations of the conventional Radon transform. 
For example, it was shown that the Radon transform on co-planar circles whose centers 
are restricted to a bounded domain is invertible \T] . Radon transforms on other smooth 
curves have also been considered [2H1]. Recently, a series of papers have explored a 
circular-arc transform which arises when the signal is generated by first-order Compton 
scattering of X-rays [5l|6]. In the imaging modality proposed in references [5l|6], the 
contrast mechanism is related to the spatially-varying efficiency of Compton scattering 
while attenuation of the scattered rays by the medium is neglected. Therefore, the 
absorption coefficient can not be recovered using this modality. 

We have recently proposed an approach which is also based on the single-scattering 
approximation but allows one to take into account and to reconstruct both the 
attenuation and the scattering coefficients of the medium [71 18]. The technique is 
applicable to either Compton scattering in the case of X-ray tomography, or to elastic 
scattering in the case of optical tomography, and was termed by us as SSOT (single- 
scattering optical tomography). The mathematical underpinning of SSOT is the broken 
ray transform of the medium which is obtained by employing collimated sources and 
detectors whose illumination/detection direction vectors lie in the same plane but are 
not on axis; the former defines the slice in which the image is reconstructed. A broken 
ray consists of two straight segments connected by a vertex and, generally, resembles 
the letter "V" (the angle in "V" can be larger than vr/2). A similar transform was 
considered recently in reference [9], where the physical mechanism of image formation 
was related to Compton scattering of gamma-rays emitted by an intrinsic radioactive 
contrast agent. 

It can be seen that generalizations of the Radon transform to non-smooth curves 
are relatively new and unexplored. This is especially true for the case of SSOT which 
employs broken rays whose directions take only few (possibly, only two) discrete values. 
This is in sharp contrast to the conventional idea that multiple ray directions are 
required for stable image reconstruction. The fact that the broken-ray transform is 
invertible is, therefore, counterintuitive. Yet, in references [7l[8], image reconstruction 
was demonstrated using a purely numerical algorithm which involved discretization of 
the integral transform and seeking the pseudo-inverse solution to the resulting system 
of algebraic equations. From the singular-value analysis of the discretized integral 
transform operator, it was found that the broken ray transform is mildly ill-posed even 
in the case when only two ray directions are used. Moreover, we have shown that, if more 
than one broken ray is used for detection, the scattering and the absorption coefficients 
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Figure 1. (Color online) Geometry a broken ray for normal incidence. 



of the medium can be reconstructed simultaneously. Another feature which is potentially 
useful is that inversion of the broken ray transform does not require multiple projections 
and can be, in principle, performed in the backscattering geometry if transillumination 
data are not available. 

In this paper, we derive and test numerically explicit image reconstruction formulas 
for the broken ray transform. In section [21 the transform and some relevant geometrical 
quantities are introduced and defined. In section |3l the image reconstruction formulas 
are derived for the cases when the inhomogeneities are purely absorptive (section 13.11) 
and when both absorbing and scattering inhomogeneities are present simultaneously 
(section l3.2p . Numerical examples are given in section HI Finally, section [5] contains a 
discussion and a brief summary of obtained results. 



2. The broken ray transform 

The physics of broken ray formation has been discussed in reference [7]. Here we focus 
on image reconstruction. The geometry of a broken ray is illustrated in figure [H for the 
particular case of trans- illumination and normal incidence. 

Imaging in SSOT is performed slice-by-slice. A slice is determined by the direction 
vectors of sources and detectors which are assumed to be sharply collimated and lie in 
the same plane. In the rectangular reference frame of figure [H the slices correspond to 
different planes x = const. Once a slice is selected, the direction vectors are fixed, but 
the positions of the sources and detectors can be scanned along the Y-axis, subject to 
the constraint that the vertex R (the ray turning point) lies within the slab < z < L. 
The position of the source is denoted by j/i and the position of the detector by 2/2, and 
every distinct source-detector pair corresponds to a unique broken ray. The measured 
intensity of a broken ray can be used to construct a data function 0(yi, ^2) as is described 
in reference [7]. The attenuation and the scattering coefficients of the medium, /i^ and 
fig, are related to the data function by the integral transform 



BR(i;2,2/i) 



I2t{rii))di-\n 



^s(R(Z/1,1/2)) 



l^s 



0(1/2, Z/i) 



(1) 
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Here fis is the background (average) scattering coefficient of the medium, the integral is 
evaluated along the broken ray BR(y2 5 Vi), r is a two-dimensional vector of position in the 
YZ-plane, and R(|/i, ^2) is the ray turning point. The problem of image reconstruction 
is to recover /Ui(r) and /is(r) from a set of data points (j){y2,yi), whereas the absorption 
coefficient fiai^i^) can be determined as the difference between /ii(r) and /Xs(r). 

If the scattering coefficient is spatially uniform and only absorptive inhomogeneities 
are to be reconstructed, it is sufficient to use a family of broken rays in which y2> Hi- 
We denote the transverse source-detector separation by A (A = ^2 — Vi)- If the angle 
between the illumination and detection directions is 6', then A varies in the interval 
(0, Amax), where Amax = -^tan^. The broken ray shown in figure [!]( a) consists of two 
segments of length Li and L2, respectively, where 

Li(A) = L-Acote = L(l-A/A^ax) , (2a) 

L2(A) = A CSC = (A/A^ax) VL2 + Al.^^ . (2b) 

If the absorption and the scattering coefficients can both fluctuate, the 
measurements of the type schematically shown in figure [T](a) are insufficient to 
reconstruct the two coefficients simultaneously. However, if the two broken rays with 
a common vortex which are shown in figure [U^b) are registered, the problem becomes 
well-determined and the two coefficients can be robustly recovered. 

The measurement scheme described above is suitable for reconstructing an image in 
the rectangular area defined by the inequalities Q < z < L and ?/min < y < Z/max, where 
L is the depth of the sample, and ymm, I/max are determined by the window in which 
the sources and detectors are scanned. In the case of optical imaging, the measurement 
scheme described above can be experimentally realized with the use of a CCD camera 
whose optical axis is tilted with respect to the normal to the slab surface. 

3. Image reconstruction formulas 

3.1. Spatially-uniform scattering 

If the scattering coefficient is spatially- uniform, /is(R) = fis, and the logarithmic term 
in the left hand side of equation ([1]) vanishes. The image reconstruction problem is then 
reduced to recovering the function fit{y, z) from the integral equation 



^JitW),z{l)]M = ct>{y^,y2) . (3) 

'BR(j/i,y2) 

It is convenient to introduce a change of variables, specifically, 

2/1 = w , 2/2 = w + A , (4) 

and define a new data function according to 
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'i/j{w,A) = (j){w,w + A) . (5) 

Thus, we parameterize the data by the position of the source, w, and by the transverse 
source-detector separation, A. The shape of a broken ray (within the slab) depends 
only on A but not on w. Therefore, we can write 

yii)=w + r]iA,i), z(£)=C(A,£), (6) 

where the functions 77(A, i), ({A, i) are independent of w. Taking a Fourier transform of 
equation ([3]) with respect to w, we obtain a generalization of the Fourier-slice theorem: 

/.Li(A)+L2(A) 

j e-^''^^^''^fit(^k,aA,i)ye = ij{k,A) , (7) 

where the Fourier transforms are defined as 

^{k,A)= ij{w,A)e''^dw, ~^t{k,z)= i^t{y,z)e"'ydy . (8) 



Thus the two-dimensional integral equation ([3]) has been reduced to the one-dimensional 
integral equation ([71) which is parameterized by the Fourier variable k. 

The integral equation ([7]) can be inverted analytically. To this end, we must specify 
the functions ri{A,i) and ({A,i). For the specific geometry of the broken rays shown 
in figure [T]( a). 



.(A£) = <!°' ^<^^(^) (9a) 

'^ ' ^ ^ ^£-Li(A)]sin^, Li(A)<£<Li(A) + L2(A) ' ^ ' 

C(A £) = <;'' ^<^^(^) (9b) 

^^ ' ^ ^ L^{A) + [i-L,iA)]cos9 , L^{A) < I < L^{A) + L^i^A) ' ^ > 

Upon substitution of these expressions into (|7j), we obtain 

rLi{A) jfcLi(A)tan6» rL 

/ fit{k,i)di+ / fitik,i)e-"'''^'''di = ^{k,A) . (10) 

Jo COS^ Jl,(A) 

We can now use the degree of freedom associated with the variable A to invert fITOl) for 
any fixed value of k. To simphfy the derivations, we introduce several new notations: 



q = kta.n6 , (11a) 

A = cot(0/2), c = cos0, /€ = — ^ = cot(^/2)cot^ , (lib) 

1 — c 

f{z) = jl{q cot e, z) , F{z) =i>(q cot 9, {L - z) tan e\ . (lie) 

The dependence of f{z) and F{z) on q is implied. Then fITU]) takes the form 
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f{i)di+-e'''' 
c 



e~""f{i)di = F{z) , 0<z<L. 



(12) 



In this equation, F{z) is known and f{z) must be found. To solve f lT^ for for a fixed 
value of q, we differentiate once with respect to z and obtain the following equation: 



-f{z) + ^e'^^ 



-iqe 



f{e)di = F'{z) , 



(13) 



K c 

where prime denotes differentiation. We then use flT^ and flT^ to find the linear 
combination G{z) = F\z) — iqF{z) in terms of f{z). The resultant equation is 



--f{z)-iq f{i)dl = G{z). (14) 

'^ Jo 

Differentiating one more time with respect to z, we obtain the differential equation 



f\z) + iKqf{z) = -KG\z) 



which has the solution 



f{z) 



-inqz 



/(O) -K e'^^^G'{i)di 



(15) 



(16) 



We then set 2; = in flT^ and find that /(O) = —k,G{0). Substituting this result into 
( 1T6l) and integrating once by parts, we arrive at the inverse solution to equation [121 



fiz) 



-K 



G{z) — inqe~ 



■iKqz 



'^''^G{i)d£ 



(17) 



One important comment on the obtained solution is necessary. The function F{z) 
in ( IT2I) is not arbitrary but such that 



F{L) 



-iKqL 



cF(0) + inq / e''"'^F{i)di 



:ii 



This can be verified directly. However, experimental measurements may result in a 
function F{z) that does not satisfy this condition. On the other hand, the inverse 
solution (TTTj) is invariant if we add to F{z) a function of the form aexp{iqz), where a 
is an arbitrary constant. It can be easily shown that any experimental function F{z) 
can be uniquely written in the form F{z) = -Frcg(^) + aexp{iqz), where F^eg{z) satisfies 
the condition ( IT8l) . Thus, inverse formula ( TT71) involves regularization, or filtering of the 
input data. 

Restoring the original notations, we find the inverse solution to ( ITOl) : 



Ht{k,z) = X 



H{k, z) — ikXe 



'ikXz 



JXki 



H{k,i)di 



(19) 



where 
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d 



H{k,z)^ [^ + tk\^p{k,A) 



(20) 



A=(L-z)tane 

The real space solution is obtained by applying the inverse Fourier transform, namely, 



Ky,z) 






(21) 



which yields 



fxtiy,z) = X 



9 ,, ,d' 



d 

ij{y, A) + K—i){y + \z, A„ 



K,{1 + k] 






^(y + K{£-A),e)di 



(22) 



A=(L-z)tant 



The above equation is the generalization of the conventional filtered backprojection 
formula to the case of the broken ray transform ([3]) . 

3.2. Spatially-nonuniform scattering 

Simultaneous reconstruction of scattering and absorption can be realized without 
making use of the entire parameter space which is available in SSOT. It is sufficient, for 
example, to use normal incidence and two different angles of detection. A particular case 
of this measurement scheme is illustrated in figure [T]^b), where the source is scanned 
along the F-axis and the intensity of two distinct broken rays is measured for every 
position of the source. The two broken rays have a common vortex at the point R and 
are denoted as BR^ and BR^. The data functions obtained by measuring the intensities 
of respective rays are denoted by (pa{yi,y2) and (pb{yi,y2)- 

At the first step, we eliminate the logarithmic term in equation [1] by taking a 
differential measurement. The differential data function is defined as 



(f>d{yi, y2) = (f>a{yi, 1/2) - Mvi^ 2/2) 



(23) 



It is easy to see that the (pd{yi,y2) can be used to reconstruct fit{y,z) even when the 
scattering coefficient of the medium is spatially-nonuniform. Indeed, because the a- and 
b-type broken rays have the same vortex, we have 



fit{yii),zii))di- / fit{yii),zii))di = Myi,y2) , (24) 

As above, the inversion formula for ( l24l) can be derived by employing the Fourier 
transform. We make the change of variables @ and define the new data function as 
ipd{w, A) = (j)d{w, w + A). The Fourier slice theorem then takes the form 



-2i 



Li(A)+L2(A) 



Li(A) 



sin kr]{A, i)]fit{ fc, C(A, i)]di = Mk, A) 



(25) 
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where ri{A,£) and ({A,i) are given by (jH])- Using these expressions, and making an 
appropriate change of the integration variable, f l25|) can be re- written as 



2i 
cose./i^(A) 



sm(k{e-Li{A))ta.ne)jlt{k,£)de = Mk,^) ■ (26) 



To solve for fit{k,z), we use the change of variables z = Li(A), (equivalently, A 
(L — z) tan 9) and differentiate twice with respect to z. This yields 



sin 6' /id 



2 



MK ^) = ^^ 1 -^^A2 + '^ ' '^''^^' ^^ 



(27) 

A=(L-2)tan6l 

The real-space solution fit{y,z) is obtained by the inverse Fourier transform. The final 
result is 



sin6' ( d'^ f°° d 

f^t{y, z) = -^ I -^ \ sgn{y - w)'il:d{w, A)dw - 2—ipd{y, A) 



■ (2^ 

A=(L-2)tane 



In the above equation, sgn(x) denotes the sign of x and it is assumed that the data 
function ip^iw^A) vanishes for \w\ > Wmax > 0, so that the integral in the right-hand 
side of (l28l) converges. 

With fit{y,z) thus determined, fis{y,z) can be obtained directly from ([T]), where 
one of the two rays (either the a-type or the b-type) is used. Finally, the absorption 
coefficient in obtained as fia{y, z) = fit{y, z) — fis{y, z). 

4. Numerical Simulations 

In this section, we illustrate the inversion formulas derived above with numerical 
examples. In the case of a constant scattering coefficient, we will use the Fourier-space 
formulas fll9j) . fl20j) and then the inverse Fourier transform (121 j) to reconstruct the total 
attenuation coefficient /Xi(y, z). In the case of spatially- varying /i^, the real-space formula 
( 128|) will be used to reconstruct Utiy^ z), then the result will be substituted into ( 12^ to 
reconstruct fis{y,z); the absorption coefficient will be obtained as fia = fJ^t — l^s- In all 
cases, the reconstructions are carried out in a rectangular area 0<2;<L, 0<2/<3L 
and the detection angle is 6* = 7r/4 so that Amax = L. 

In what follows, the optical coefficients are decomposed as 

lit{y,z) = Jit + 5nt{y,z), Hs{y,z) = iJ,s + 5ns{y,z), Haiy , z) = Jia + 5 Ha{y , z) . (29) 

Here the quantities with a bar are constant background values of the respective 
coefficients and the (5-terms represent the inhomogeneities. In performing image 
reconstruction, it will be assumed that the support of the functions 6fia{y,z) 
and 6ns{y,z) is completely contained in the rectangular area where the image is 
reconstructed. This assumption is, strictly speaking, violated if Gaussian targets are 
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used, as will be done below. However, the Gaussian inhomogeneities are exponentially 
small outside of the imaging area, and the error incurred due to the inaccurate 
assumption is also exponentially small. 

4.I. Spatially-uniform scattering 

We consider first the case of a medium with Sfig = and a spatially-inhomogeneous 
attenuation coefficient (due to spatially-varying absorption) in the shape of a square 
with sharp boundaries, viz, 

Ml/,^) = |^*' l^-^o|<a/2and|.-.o|<a/2 _ 
^^^' ' [ , otherwise ^ ' 

Thus, the attenuation coefficient inside the square is twice the background value of 
fit- The square is centered so that yo = L, zq = L/2 and its side length is a = L/2. 
The Fourier-space data function ilj{k,A) was obtained analytically from ( TT9l) by direct 
integration. Then the variables A and k were sampled and Ht{y, z) was reconstructed 
using discrete samples of ?/'(A;,A). More specifically, the samples have been used 
to evaluate numerically the equations ( 12T|) and flT9|) . The variable A was sampled 
as A„ = hn, where n = 0,1, ... N and h = A^ns^^/N. Note that for 6 = ir/A, 
Amax = L. The Fourier variable k was sampled in the interval [—ir/h, ir/h] according to 
km = (ti /h)(m/N — 1), where m = 0, 1, . . . , 2N . The derivative in (120|) was computed 
as the central difference 



dti>{k,A) 



DA 



i}{k,An+i)-i'{ k,An-i) , 

, n = 0,l,...,N . (31) 



2h 

A=A„ 



In applying the above formula, the boundary condition ^{k,A_i) = ■ip{k,A]\r^i) has 
been used. Note that the substitution A ^ {L — z) tan 6 contained in the formula 
(1201) did not require any re-sampling because, in the geometry used, tan^ = 1. More 
generally, however, re-sampling and interpolation is required to apply f l20|) to a uniformly 
sampled function ipi^, A). The integral over i in ( TT9l) was evaluated numerically using 
the trapezoidal rule. Finally, the image was reconstructed on a rectangular grid with 
the same step as was used to sample the variable A, that is, h. Reconstructions of the 
total attenuation coefficient ^t{y, z) obtained as described above is shown in figure |2] for 
two different values of the parameter N . It can be seen that the reconstruction contains 
artifacts. When the number of samples, A^, is increased by the factor of 10, the support of 
the artifacts is reduced by the same factor, but the amplitude stays unchanged. We have 
verified that the L^ norm of the discrepancy ^{y, z) = ^i {y, z) — ^^^'^'^°^^ '^^^ ® ' (^y^ z) 
tends to zero when A^ — ?■ oo, yet the maximum amplitude of the relative error, 
max[^'^*™'')(?/, z)/ iJit{y, z)], remains of the order of unity. 

The relatively low image quality seen in figure [2] is caused by the sharp discontinuity 
of the target. This statement is confirmed by considering a smooth target in which the 
attenuation coefficient is given by a Gaussian, viz. 
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Figure 2. Reconstructed dfit{y,z) for the square target (|30p using different number 
of samples A^, as labeled. 



Sfit{y,z) = fitexp 



{y - Vof + {z- z^f 



a'' 



(32) 



The data function for this target, tlj{k,A), can be computed analytically (the result 
contains the error function). We have used the same sampling procedure as above and 
a discretized image of the target (132|) has been reconstructed using equations ( 120|) . (fT9|) 
and (!2T]l . The results are shown in figures [3] and HI It can be seen that an accurate 
quantitative reconstruction is obtained for the smooth target of the type (1321) . When the 
linear scale is used to represent the data, no artifacts are visible in the reconstructions 
and the reconstructed data coincide quantitatively with the model. However, when the 
logarithmic scale is used on the vertical axis, as is done in panels (b,d) of figure HJ 
the artifacts become clearly visible. The amplitude of the artifacts in the lateral cross 
section of the image (figure |4](b)) is, in fact, less than 2% of the plot maximum. 



4-2. Spatially-nonuniform scattering 

We next consider the case when both the scattering and the absorption coefficients of 
the medium are varying. We have modeled the inhomogeneities as Gaussians 



Sfis,a{y,z) =/is,aexp 



{y - y3,af + iz- Zs^af 



a^ 



(33) 



where {ys, Zs) and {ya, Za) are the centers of the scattering and the absorbing 
inhomogeneities. Note that the amplitude of each coefficient in the center of an 
inhomogeneity is twice the background value. Both kinds of inhomogeneity were 
assumed to be present in the medium simultaneously but not overlap. Thus, the 
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a = 9h (T = 15h a = 21h a = 30h 

Figure 3. Reconstructed 5iit{y,z) for the Gaussian target ([5^ using N = 120 and 
different values of a, as labeled. 



scattering inhomogeneity was centered at the point {ys = 3.125L,Zs = 0.5L) and the 
absorbing inhomogeneity was centers at (y^ = 0.875L,Za = 0.5L). We have considered 
different values of the parameter a. However, in every reconstruction, a was the same 
for the absorbing and the scattering inhomogeneity. 

In this subsection, we have used the real-space image reconstruction formula f l28|) . 
The data function was obtained by analytical integration and the image reconstruction 
was performed by sampling the variables w and A in the data function ipdiWjA) on a 
rectangular grid with the step h, as is described in more detail in the previous subsection. 
In the reconstructions of this subsection, A^ = 120 and h = L/N. The derivatives in 
fl28|) were computed by central differences and the integral by the trapezoidal rule. 

We consider below two cases. The first case corresponds to /i^L = 2.4 and 
jlaL = 0.24, so that the background scattering coefficient is ten times smaller than the 
background scattering coefficient. The challenge here is to reconstruct the absorbing 
inhomogeneities in the presence of much stronger scattering inhomogeneities. In the 
second case, /x^L = fiaL = 2.4, so that the strength of the absorbing and the scattering 
inhomogeneities is the same. Image reconstruction for the first case is illustrated in 
figures [5] and [61 Here we plot total optical coefficients (including the background) rather 
than the fluctuating parts 6 fit, etc., as was done in figures [2]l3]l4l It can be seen that 
a good image quality and a quantitative agreement with the model are obtained for 
the scattering and the attenuation coefficients. The reconstruction of the absorption 
coefficient is not as good. This is because the relatively small quantity fia was obtained 
by finding the numerical difference between the two much larger quantities fit and /x^. In 
the case when the magnitudes of the scattering and the absorbing inhomogeneities are 
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Figure 4. (Color online) The data of figure [3] are shown here as cross sections along 
the straight lines drawn through the center of the inhoniogeneity in the directions 
parallel to the F-axis (a,b) and to the Z-axis (c,d). The panels (a,c) use the linear 
scale and the panels (b,d) use the logarithmic scale for the vertical axis. Dots represent 
the reconstructed values and the model function is represented by the dashed lines. 
Every second reconstructed data point is shown in semi-logarithmic plots. Some of the 
data points are not shown in the semi-logarithmic plots because the respective values 
are negative or too small to be displayed. 
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Figure 5. Simultaneous reconstruction of tire absorbing and scattering coefficients 
for the case jlgL — 2.4, flaL = 0.24L. The three columns represent the attenuation, 
scattering and absorption coefficients, as labeled. The letter "M" indicates "model" 
and the letter "R" indicates reconstruction. The first two rows correspond to cr = 21h, 
the next two rows correspond to a — 9h and the last two rows correspond to cr = 3/i. 
Here h = L/N and N = 120. Every plot is normalized to its own maximum. 



the same, a quantitatively accurate reconstruction of all three coefficients is obtained, 
as is illustrated in figures [7] and [HI 

5. Conclusion 



We have derived and tested numerically image reconstruction formulas for the broken 
ray transform (also referred to as the V-line transform). The obtained formulas 
are generalizations of the filtered backprojection formula of the conventional Radon 
transform. We have shown that the broken ray transform is in certain aspects more 
useful as it allows one to reconstruct the scattering and the absorption coefficients of 
the medium simultaneously. 

Inversion of the broken ray transform is mildly ill-posed. The ill-posedness strongly 
affects image reconstruction when the target has sharp boundaries. We conjecture that 
the ill-posedness can be regularized by using additional rays. 

The broken ray transform considered in this paper may become useful on its own 
merits in the situations when scattering in the medium is not negligible. However, it 
can also prove useful in conjunction with the conventional Radon transform in which 
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Figure 6. (Color online) The data of figure [5] for the case a = 21h {h = L/N 
and N — 120) are shown here as cross sections along the straight line z — Za = Zs 
which intersects the centers of the absorbing and the scattering inhomogeneities. Total 
attenuation and scattering are shown in panel (a) and absorption is shown in panel 
(b). Centered symbols correspond to the reconstructed values and the dashed lines 
correspond to the model. Every second reconstructed data point is shown. 
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Figure 7. Same as for Fig.[5]but for figL — jiaL — 2.4. 
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Figure 8. (color online) Same as for Fig.[S]but for figL = jiaL = 2.4. 

the intensity of ballistic rays is measured. Since the ballistic ray is a special case of 
the broken ray, it seems plausible that a reconstruction formula can be obtained which 
would utilize the measurements of ballistic and broken rays simultaneously. This will 
be the subject of our future work. 
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